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D , Abstract 

In this paper we consider a charged massless scalar quantum field 
operator in the spacetime of an idealized cosmic string, i.e., an in- 
^ ' finitely long, straight and static cosmic string, which presents a mag- 

■ netic field confined in a cylindrical tube of finite radius. Three distinct 

situations are taking into account in this analysis: i) a homogeneous 
field inside the tube, ii) a magnetic field proportional to 1/r and Hi) 
a cylindrical shell with ^-function. In these three cases the axis of the 
infinitely long tube of radius R coincides with the cosmic string. In 
order to study the vacuum polarization effects outside the tube, we 
explicitly calculate the Euclidean Green function associated with this 
system for the three above situations, considering points in the region 
outside the tube. Having these Green functions we calculate the renor- 
malized vacuum expectation values, (4>*(x)&(x))n en . and (T^(x))# en ., 

*E-mail: spinelly@fisica.ufpb.br 
^E-mail: emello@fisica.ufpb.br 



1 



associated with the charged field. In the evaluation of these vacuum 
polarization effects, two contributions appear for the three models. 
The firsts are the standard ones due to the conical geometry of the 
spacetime and the magnetic flux. The seconds contributions appear as 
an extra term. They are corrections due to the finite thickness of the 
radius of the tube. These extra terms provide relevant contributions, 
even for points very far way from the system, like a long-range effect. 
PACS numbers: 98.80.Cq, ll.10.Gh, 11.27.+d 
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1. Introduction 



One of the firsts models which present topological defects was investigated 
by Nielsen and Olesen [1] several years ago. In their analysis topologically 
stable solutions named vortices, appears as a consequence of spontaneous 
breakdown of local Abelian and non-Abelian gauge symmetries in a classical 
field theory. These solutions correspond to infinitely long objects which carry 
magnetic fluxes. Later, Garfinkle investigated the Abelian model considering 
its influence on the geometry of the spacetime [2]. He found, as in flat 
spacetime, static cylindrically symmetric solutions representing vortices. He 
also shown that asymptotically the spacetime around the cosmic string is a 
Minkowiski one minus a wedge. Its core has a non-zero thickness, and the 
magnetic field vanishes outside of it. Two years later Linet [3] obtained, under 
some specific condition, exact solution for the metric tensor. He was able to 
show that the structure of the respective spacetime corresponds to a conical 
one, with the conicity parameter being expressed in terms of the energy 
per unity length of the vortex. The complete analysis about the behavior 
of the gauge and matter fields near the cosmic string's core can only be 
obtained numerically. Some recent numerical analysis [4] about the structure 
of supermassive cosmic strings show that two different kind of solutions for 
the metric tensor exist. 

A complete analysis about the behavior of a charged field in the neigh- 
borhood of a U(l)— gauge cosmic string, must take into account not only 
the influence of the geometry of the spacetime but, also, the influence of the 
magnetic field. Two distinct cases can be analysed: i) The first one is to con- 
sider the string as an idealized infinitely thin linear topological defect, having 
a magnetic field running along it. This case can be treated analytically, ii) 
The second approach is to consider the non-zero thickness for the string. 
Analytically this problem becomes intractable. An intermediate approach is 
what we shall adopt here. We shall use an approximated model consider- 
ing the spacetime produced by the string as being conical everywhere, but 
having a magnetic field inside a tube of non-zero thickness surrounding it. 
In this way some improvement is introduced when compared with the ideal 
case. Moreover, we shall admit specifics spatial behaviors for the magnetic 
field inside the tube trying to approach the real situation. 

The vacuum polarization effects due to a magnetic field confined in a 
tube of finite radius in Minkowski spacetime has been first analysed by Sere- 
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bryanyi [5]. A few years later this effect in a charged massless scalar field on 
a idealized cosmic string spacetime, has been calculated by Guimaraes and 
Linet [6]. There a magnetic flux running through the line singularity was 
considered. As a consequence the renormalized vacuum expectation value 
associated with the energy- momentum tensor, (T^(x)) R en ., presents a con- 
tribution coming from the geometry of the spacetime and the magnetic flux. 
In order to obtain this result the respective Green function was determined. 
In a particular choice of gauge, this Green function depends explicitly on 
the factor 7, the fractional part of the ratio of the magnetic flux <fi by the 
quantum flux, 0o — 27r/e. 

More recently Sriramkumar [7] has calculated the vacuum fluctuation of 
current and energy densities for a massless charged scalar field around an 
idealized cosmic string carrying a magnetic flux. His approach to obtain the 
Green function takes into account the presence of the vector potential in 
the differential operator, = <9 M — ieA^, which presents the advantage to 
calculate the two-points function without imposing any additional boundary 
condition on the field. Also the vacuum polarization associated with massless 
scalar field in a infinitely thin cosmic string embedded in a Randall-Sundrum 
type 2 brane world, has been calculated [8] . In lowest order of perturbation, 
{T^Ren., consists of two parts. The first one depends only on the conicity 
parameter and the second is a correction which depends on the parameter k 
connected with the cosmological constant. 

Allen, Kay and Ottewill [9] have analyzed the vacuum polarization effect 
of a massless scalar field on a realistic cosmic string spacetime, considering 
generically the effect of the string's core through the non-minimal coupling 
between the scalar field with the geometry: £1Z, being TZ the Ricci scalar 
and £ the coupling constant 1 . In this case the respective Green function 
presents two parts: The first one, named "regular" by the authors, coincides 
with the usual one obtained in the literature for an idealized cosmic string, 
see [14] for example. The second contribution comes from the nonzero core's 
radius of the string. Although the radius of the cosmic string's core is very 
small 2 , the effect of this finite thickness on the vacuum polarization of a 
massless scalar field, is a long-range effects. Moreover it is of the same order 

1 In [9] was not considered a specific model for the inner structure of the string; however, 
it was considered be infinitely long, straight and static. 

2 For cosmic string formed by a typical grand unified theory is of order 10~ 30 cm. 
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of magnitude as contribution coming from the geometry up to distance from 
the cosmic string that exceed the radius of the observable Universe. 

Motivated by this result we raised a question about the consequence on 
the vacuum polarization of a charged massless scalar filed in the U(l)— gauge 
cosmic string spacetime considering this string surrounded by a magnetic filed 
confined in a tube of finite radius. As in the Allen et al paper, we expect 
some consequence of the finite thickness of the tube on the vacuum polar- 
ization effect. In order to investigate this subject we decided in this paper, 
for simplicity, to analyse the behavior of a quantum charged massless scalar 
field in a idealized cosmic string spacetime, considering the presence of a 
infinitely long tube of finite radius R of magnetic field. Three distinct cases 
were taking into account: 

i) a homogeneous magnetic field inside the tube, 

ii) a magnetic field proportional to 1/r inside and 
in) a magnetic field on a cylindrical shell. 

In all these three cases the axis of the infinitely long and straight tube coin- 
cides with the cosmic string. 

The idealized model for a infinitely long straight static cosmic string 
spacetime can be given in cylindrical coordinates by the line element below: 

ds 2 = -dt 2 + dr 2 + a 2 r 2 d6 2 + dz 2 , (1) 

where q is a parameter smaller than unity which codify the presence of 
a conical two-surface (r,0). In fact for a typical Grand Unified Theory, 
a = 1 — O(10 -6 ). (Here we considered that this cosmic string is lying along 
the z-direction.) 

A magnetic field along the z-direction can be represented by H(r) = 
H(r)z. Now assuming that the field has a finite range in the radial coordinate, 
we are particularly interested in the three models below: 



i) H(r) = — — r), homogeneous field inside (2) 

ii) H{r) = — - — Q(R — r), field proportional to 1/r inside (3) 
Hi) H{r) = — ^— S(r — R), cylindrical shell, (4) 
where R is the radial extent of the tube, and <fi is the total flux. The ratio of 
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the flux by the quantum flux O , can be expressed by 0/0o — N + 7, where 
N is the integer part and < 7 < 1. 3 

Considering both effects, due to the geometry of the spacetime and the 
finite thickness of the tube of magnetic field on the vacuum polarization of 
a charged massless scalar field, we observe that there appear two distinct 
contributions for each case: For the cases (i) and (ii), the firsts ones are 
due to the geometry and the magnetic flux. They are similar to the result 
found by Guimaraes and Linet [6]. The seconds contributions are due the 
nonzero radius of the tube, R. The latter present similar dependence with the 
radial coordinate. In fact their dependence on the ratio R/r is of power type 
besides an overall 1/r 4 dependence. These results indicate that the detailed 
information about the structure of the magnetic field inside the tube is not 
so relevant. The third case is completely different from the others. The 
second contribution to the vacuum polarization is now very similar to the 
result obtained in [9]. It presents a logarithmic dependence on the radial 
coordinate, being a long-range contribution. 

Taking into account that the topological defects are of evident physical 
interest (they are considered as good candidates to explain some components 
of anisotropy on the cosmic microwave background [10], gamma ray burst 
[11], gravitational waves [12] and highest energy cosmic ray [13]) it seems 
reasonable to investigate the above mentioned problem in more detail. 

This paper is organized as follows. In the section 2, we explicitly calculate 
the Green functions associated with this system for the three different mod- 
els. Having now these Green functions, in section 3 we explicitly calculate 
the renormalized vacuum expectation values, (&*(x)&(x)) R en ., in complete 
form. Using a judicious criteria of approximation we present the behavior for 
these quantities. In section 4, we calculate the renormalized vacuum expec- 
tation values of the energy-momentum operator, (T"(x))R en .. Using the same 
procedure as in the previous section we present a approximated expressions 
for these results. We leave for section 5 our conclusions and most relevant 
remarks about this paper. Finally we want to say that we shall use in this 
paper h = G = c = 1. 

3 These three models for the magnetic field has been considered previously by Bordag 
and Voropaev [15] to analyze the quantum mechanical behavior of spin— 1/2 and gyro- 
magnetic g ^2 particle in a field of a magnetic string. 
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2. Green Function 



The Green function associated with a charged massless scalar field must obey 
the non- homogeneous second order differential equation below 



-.V, 



-gg^V u ] G(x,x) = -5 {A \x - x'), 



(5) 



where = 

Because it is easier to calculate the vacuum average of the square of the 
field and the energy-momentum tensor operators using the Euclidean version 
of the Green function, let us solve the above equation adopting the Euclidean 
extension of the cosmic string metric tensor given in (1), by performing a 
Wick rotation t — > it. Now if we want to represent a magnetic field along 
the z-direction, the four-vector potential reads 



A„ = (0,0,A(r),0) . 



(6) 



So the equation (5) becomes 
d 2 1 d ( d \ d 2 1 



dr 2 r dr \ dr 



dz 2 



a 2 r 2 



dd 2 



G(x,x') = 



-—5(r - r')5(r - r')5(9 - 6')5(z - z') . (7) 



ar 



Now, in order to reproduce the magnetic fields configurations given by 
(2)- (4), we have only to write the third component of the vector potential by 



2txol 



a(r) 



(8) 



For the two firsts models, we can represent the radial function a(r) by: 
a(r) = f(r)Q(R - r) + 0(r - R) , (9) 

with 

r 2 /R 2 , for the model (?) and 



m = 

For the third model, 



r/R, for the model (ii). 
a(r) = e(R-r). 



(10) 
(11) 
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Due to the cylindrical symmetry of this system the Euclidean Green function 
can be expressed by 

I l ^~ > r oo r oo 

G(x,x') = £ e m ^ / dke ik ^ / dwe™^ g n {r,r>) . (12) 

W 6 n=-oo J -™ J ~°° 

Before to specialize on the specific model let us write down the non- 
homogeneous differential equation obeyed by the unknown function g n (r, r'). 
Substituting (12) into (7) and using the standard representation to the delta 
functions in the temporal, angular and ^-coordinates, we arrive to following 
differential equation for the function g n (r,r'): 



+ -4- ~ - 2neA + e ^ 2 ) - 1 

dr z r ar a z r z 



g n (r,r') = — —5{r - r') , 
ar 



(13) 

where f3 2 = k 2 + u 2 . 

So, all the informations about the magnetic field, the radius of the cylin- 
drical tube and the conicity of the spacetime are contained in g n (r, r'). How- 
ever, before to obtain specific solution, corresponding to specific model of 
magnetic field, some general properties must be satisfied by this function. Let 
us define by g<(r,r') the solution of (13) regular at r — > 0, and by g>(r,r') 
the solution that vanishes at infinity. These two solutions must satisfy the 
continuity condition at r = r', with their first derivative discontinuous at 
this point. 

It is of our main interest to investigate the vacuum polarization effect for 
external points to the magnetic flux. As we shall see, because the information 
about the field in this region must taking into account, the above mentioned 
effect will present contribution coming from the inner region. Now, let us 
consider solutions of (13) with both r and r' greater than R. 

Let us first consider the models (i) and (ii). The inner solution, g<(r,r'), 
of (13), corresponds to r < r'. However two regions must be considered 
separately: from to it! and from it! to r . Integrating out in region r < R, 
we have: 

< ( a _ / A i H i( r ), for r < i? 

9n[r ' r) -\ BiJ^P^+C^iPr), foTR<r<r', 

where 

(n — 5) , , 

y = - (15) 

a 
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being 

S=^ = N + ,. (16) 

Hi(r), for % = 1,2, represents the solution associated with the two firsts 
models: 

*.<-> = \ M nM (^) , (17, 
with (n = (2 - ^2)/2 and Ai = n/2«, and 

if 2 (r) = i=M (T2 , A2 (Cr), (18) 
V r 

with <r 2 = ^(5 2 - (3 2 R 2 a 2 )- 1/2 , A 2 = n/a and C = ^(5 2 + plPa 2 ) 1 ' 2 . 

In both definitions above, Eq.s (17) and (18), M a> \, represents the Whit- 
taker function [16]. The continuity condition of g„ and its first derivative at 
r = R, allow us to obtain the ratio D l n = C l n / ' B l n as a function of (5R only: 

H'(R)IU(5R) - HAR)IU(3R) 
n[P ) H^K'^PR) - H>(R)K H ((3R) ■ 1 ] 

In all the expressions, l\ v \ and K\ v \ are the modified Bessel functions. 
The outer solution of (13) is given by 

g>{r,r') = W n K\ v ^r\ for r > r. (20) 

All the above constants are determined by imposing the boundary con- 
ditions explained before on g< and g> at r = r'. In this way they will be 
function of r'. Finally we obtain the result to the complete radial Green 
function below: 

9n(r,r') = -LK W \(p r> ) (/ M (/3r<) + D\K W \(p r< )) for r and r > R. (21) 

In the above equation r> (r<) is the bigger (smaller) value between r and r'. 
Substituting this solution into (12) and after some intermediate steps, the 
Euclidean Green function acquires the following expression: 

e iNA0 e iAe sinh(7M /«) + sinh[(l -7)w /a] 
87r 2 arr' sinh u cosh(-u /o;) — cos A6 



+ 4tt 2 



oo 

]T e- Ae ^(^)K| H (/5r)ir| H (/3r'), (22) 



n=— oo 



where 

r 2 + r ,2 + (Ar) 2 + (A^) 2 

coshw = — ^ . (23) 

2rr 

We can observe that the first term of the above Green function is, up 
to a gauge transformation, equivalent to the result presented by Guimaraes 
and Linet [6] for the massless scalar field. The second term presents, besides 
the dependence on the a parameter and the magnetic flux, a dependence 
on the radius R of the tube through the constant D l n . Moreover, this Green 
function obeys the boundary condition G(0' = 0+2tt) — G(0' — 0), which is a 
consequence of having used explicitly the vector potential in the differential 
equation (5). Although we have calculated this Green function for points 
outside the magnetic tube, there remains contributions coming from this 
region that do not disappear by taking 7 = and A0 = 2n. 

As to the third model, the solution to g n (r,r') is: 

a<(rr') = l A^ M (/H for r < i? 

9n {r ' n \ B n I\ nH (Pr) + C n K lnH ((3r), for R < r < r , ^ 

and 

g>(r,r') = E n K\ n/a \{f3r), for r > r . (25) 

Once more, the complete solution to this radial function can be obtained 
adopting a similar procedure as the previous case. This solution is: 



g n (r,r) = --iif| n |/ a (/3r>) (i|„|/ a (/3r<) + D n K\ n \ /a ((3r < 

for r and r > R. (26) 



Again imposing the continuity condition on (24) and its first derivative at 
r = R we get D n = C n j B n as function of (3R: 
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Finally adopting similar procedure as in the previous case, we obtain: 



G(x,x') 



87v 2 arr' sinh-u 



sinh(-u / a ) 



cosh(w /a) — cos A6 



+ ^ jT d/3/3J (/?V(Ar) 2 + (Az)- 



x 



£ e mM D n ((3R)K H/a ((3r)K H/a {{3r') 

n=— oo 



(28) 



For this case we can see that the first term in this Green function coincides 
with results found by Smith [14] and by Linet in [17] for a massless scalar 
field. Moreover, the second term contains the information about the radius 
of the tube through the constant D n . 

So having the Green functions our next steps will be the calculations of 
the vacuum expectation values associated with the modulus of the scalar field 
and its energy-momentum tensor. 



3. Computation of {&(x)<h(x)) Ren . 

In this section we want to calculate the renormalized vacuum expectation 
values associated with the modulus square of the scalar field in this physical 
system for the three distinct configurations of magnetic fields. As we shall see 
these vacuum average will depend on the conicity parameter a, the fractional 
part of 0/0o and the radius of the tube, R. Because the two firsts models 
present similar Green functions, we shall develop initially the calculations 
considering these two models together and later the third one. 

The renormalized vacuum expectation value of the modulus square of the 
charged massless scalar field on idealized cosmic string spacetime is given by 

(®*(x)$(x)) Ren . = lim [G(x, x') - G E (x, x')] , (29) 

x —*x 

where Ge(x,x') = 4^27~7p ; is the Euclidean Green function for a flat four- 
dimensional spacetime, and G(x, x') the Euclidean Green function associated 
with this system under investigation. 

As we have explained before, we shall develop first this calculation for 
the models (i) and (ii), whose Green functions are given in (22) with D % n 
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given by (19). A very interesting aspect of these Green functions is that the 
correction due to the magnetic tube's radius is finite in the coincidence limit; 
so the result is: 



mxMx)) Ren . = ^- i -^[l-a 2 -6 7 (l- 7 )] 

1 

+ 



—J d/3/3 ]T DU^Kf^r). (30) 

JO n=-oc 



In order to provide a quantitative information about the behavior of the 
correction in above equation, we must adopt the following procedure: i) 
Redefining a new variable v — (3r, it is possible to see that the second term 
presents on overall 1/r 2 dependence, the same as presented by first term. 
However there remains an extra radial dependence on the integrand through 
the coefficient D l n (vR/r). So the complete information about the second 
term can only be obtained numerically specifying the ratio R/r. ii) An 
approximated behavior can be exhibited if we analyse the integrand of (30). 
Near the origin K\ v \(v) behave as v _2 I"I, on the other hand D l n {vR/r) present 
a dependence v 2 ^ as we shall see later. So the product vD l n (vR/r)K, 



M 



v 



vanishes linearly with v. For large v, vK 2 u ^{v) decays as e 2v . Because 
D l n (vR/r) grows slowly than e 2vR l T for large v, its most significant relevance 
in the integrand happens only up v of order unity. However because we are 
mostly interested to study the vacuum polarization for points very far from 
the cosmic string, we shall consider R/r <C 1 in the our analysis. So we may 
approximate the integrand of (30) assuming to the coefficient D l n its first 
order expansion. This is equivalent to consider the coefficients D l n (f3R) with 
f3R 1. Taking in (19) f3R 1, the most relevant term for both cases are 
given by 




where 



D '^ R) = -rOHTTjrM ' H ^ > (31) 



S" = {— - Ma„a,0V«) + (l + -) M, uXl (S/a), (32) 



a } \ a / 

f? = M XlM (5/a), (33) 
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92 = - 5) M X2M {25/a) + (I + ^) M 72iA2 (25/«) , (34) 

and 

K = M X2 , Xl (26/a) , (35) 

being 7 i = (n + 2a) /2a, 72 = (n + a)/a and z/ = "~"^~ 7 . 

Using now these approximated expressions to the coefficients D l n , and 
with the help of [16], Eq. (30) becomes: 



mx)Hx)) Ren . = ^J-^[l-^-6 7 (l- 7 )] 

1 < 





u\ 




1 + 2 


\v\ 




- (36) 



(2tt)W 2 n f^ 1 + 2\u\ \gf + \v\fpj V r 7 



The first term of the above result depends only on the fractional part 
of 0/0o, besides the parameter a. On the other hand the second term, the 
correction, depends on the integer part of the the ratio <p/4>o, N. However, 
because R/r is assumed to be smaller than unity, the most relevant correction 
comes from the n = N. In this case the extra radial dependence of the second 
term is dominated by (R/r) 2 " 1 ^. 

As we have already explained our approximation consists to substitute 
into (30), the expression (31) to the coefficients D l n . In fact doing this other 
contributions given by the complete expression (27) are missing; only the 
first order in the expansion was kept. In order to be sure that our approx- 
imation is a good one, we also analysed numerically the integrand of the 
second contribution of (30) using the exact expression (27) and the approx- 
imated result (31)-(35) for the coefficients for both cases i = 1 and 2. For 
simplicity only we take in the summation N = 0, for both cases. This ap- 
proximation is once more justified by the figures 1(a) and 1(6) which exhibit 
comparatively the integrand of (30) written in terms os the new variable 
v, I{y) = vY.D % n {yR/r)K^{y), for the models (i) and (ii), respectively, 
considering the complete expressions for the coefficients and their respective 
approximations. We can see a good agreement between both curves for each 
case. (In these numerical analysis we have considered only the term n = 
in the summation for both models. Considering other terms, for example n 
running from —20 to 20, we have not observed any significant difference when 
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compared with the n = term. So the summation in (36) can be mainly 
represented by its most relevant contribution.) 

So accepting this fact as general, our final result to the vacuum polariza- 
tion effects becomes: 



1 — a 2 — 67(1 — 7) 



487r 2 r 2 a 2 | 

12 7 
0(27 + a) 



a' 



9, 



N 



7/«/, 



N 



9? + 7/«/, 



N 




(37) 



From the above expression, it is possible see that the second term has the 
same order of magnitude as the standard one, up to ratio r/R = |x| a//27 , 
being x given by 



x 



a 



12 7 (gf-i/afi 



v- 



1 _ a 2 + 67 ( 7 _ 1) a + 2 7 + 7/a/ j 



iV 



(38) 



For smaller ratio the second contribution becomes more relevant than the first 
one. Once more, by numerical analysis we observe that for N = 0, a = 0.99 
and 7 = 0.3, there happen a cancellation between both contributions for r/R 
of order 5 for both models. For N — 1 and assuming the same values for 
a and 7, the seconds terms changes their signs and there is no cancellation 
anymore. 

After this discussion about the two firsts models, let us analyze the third 
one where the respective Green function is given by (28). For this case the 
renormalized vacuum expectation value of the modulus square of the scalar 
field becomes: 



+ 



A8n 2 r 2 a 2 
1 



a 







d(3(3 ]T D n ((3R)K 2 nl/a ((3r) , (39) 



with the coefficients D n (f3R) given by (27). Although the physical system 
that we are analysing is completely different from the one analyzed by Bruce 
Allen et al in [9] , there appears a strong similarity between our result with 
theirs. Again if we are inclined to adopt an approximate expression to the 
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coefficients D n considering (3R <C 0, two distinct behavior are obtained: 
For n^O, 

2 ( \v\-\n\/a \ (f3R\ 2lnl/a 

D n {(iR)- r (| n |/ a + i) r (| n |/ a ) [\u\ + \n\/a) [ 2 ) ' (4 ° } 

and for n = 0, and 

= S/a [to<M) + c}-l • <41) 

where C is the Euler constant. As we can see for n ^ 0, D n (/3R) vanishes 
at least as fast as (f3R) 2 ^ n ^ a when (3R — > 0. On the other hand, D Q ((3R) 
vanishes only with the inverse of the logarithm, so slowly. In this way the 
most relevant contribution to the summation in (39) comes from the n — 0. 
This same conclusion was previously obtained by Bruce Allen et al for another 
physical system. There they used a more convenient notation to express this 
coefficient, which also we adopt here: 

D ^ H) = M/b (42) 

where in our case 

q = |e^. (43) 

So the most relevant contribution to the correction of the renormalized 
vacuum expectation value of the square of modulus of the scalar field in (39) 
is proportional to: 

I 00 dBB^Mll = I r dvv KI{V) (44) 
io PP H(3/q) Wo dV \n(v)-\n(qr)- [ ^ 

As it was pointed out in [9] this approximation presents a problem because 
the above integrand has a pole at v — qr. However, this pole is a consequence 
the approximation adopted, and it occurs when the argument of D is much 
greater than unity, where the approximation is no longer valid. Moreover 
analyzing the behavior of the integrand, they observe that vKq(v) vanishes 
at v — 0, peaks around v of order 1, and decays as e~ 2v for large v. On 
the other hand D vanishes for v — and grows less slowly than e 2vR ^ r 



15 



for large v. So, finally they conclude that (44) can be well approximated 
discarding the term ln(v). (Also a numerical analysis in their paper supports 
this approximation.) 

Here in our case a similar conclusion can be reached. In fact choosing 
7 = 0.2, N = and a = 0.99, qr = ^e~ c+a/s in of order than 10 5 for 
r/R = 10 3 , so f3R ~ 10 2 . Consequently the pole occurs in the region where 
the approximation is no longer valid. We shall not repeat their argument 
to accept the same approximation procedure; however we present in Fig. 
2 a numerical evaluation of this approximation, comparing the integrand 
of the second contribution of (39) using the exact expression to its more 
relevant term, D , given by (27) and its approximated expression given in 
(44), discarding the ln(t>) term. So we can conclude, supported by those 
evidences, that the second term of (39) can be approximated by 



1 1 



/ dvvKl(v). (45) 
Jo 



(2n) 2 r 2 ln(gr) jo 

Finally with the help of [16] we integrate the above equation, providing to 
(39) the expression below: 

1 1 2 1 

<*•(*)*<*»»- = iSv^" - s^Mie-^) ' (46) 

From the above expression we can see that the explicit dependence on 
the flux tube appears only on the second term. Taking R — > or 5 — > 0, this 
term vanishes. (The first term contains only the contribution coming from 
the geometry and vanishes when a — > 1). 

As in the physical system analyzed in [9] , there is in this long-range 
effect in (46) coming from the special structure of the magnetic tube around 
the cosmic string. 

Because a« 1, (46) can be written as 



{a- 1 - 1) ln(f e" c + Q /5) 



(47) 



where we see that the second term inside the bracket vanishes slowly as 
r — > oo. In fact it can be noticed that the second term is greater than the 
unity up to 

r I r 3 ni ' 

(48) 



r 1 

- exp 



R 



a 1 — 1 o 
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which a = 0.99 and 5 = 0.2 we have r/R = 10 127 . As in [9], even for R = 
10~ 30 cm, the correction in (47) becomes more relevant than the geometrical 
effect up to r xs fO 97 cm which exceed the radius of the observable Universe. 

4. Computation of (T^(x))n en . 

The renormalized vacuum expectation value of the energy-momentum tensor 
associated with this system under investigation is given by: 

(fp Ren . = KmjD^Gix, x') - D»; (m G E (x, x')] , (49) 
with the operator D^, a ^ expressed as 

= (1 " 20^X " Z(D,>ir + D»>D") + (2£ - \)6;D a D a ' , (50) 

being D a = V CT — ieA a and the bar denoting complex conjugate. Again, in 
order to take into account the presence of the three magnetic field configura- 
tions given in the beginning of this paper we write the four vector potential 
in the form A a = (0, 0, 4^-, 0), with a(r) being given by (9) and (10), for the 
firsts two cases and by (11) for the third one. The respective Green functions 
G(x,x'), are given by (22) and (28). 

Because we are convince ourselves about the approximation procedure 
adopted in the last section, here we shall consider in all the calculations only 
the leading term of the coefficients D l n for each case. 

Let us start with the two firsts models. After some intermediate calcula- 
tions, we arrive at: 

(f°(x)) Ren . = (f°(x)) Reg . + -L- jT d(3(3Di =N {(3R) \2(3 2 - ±) x 

w + ^ 2 + 4 (e - £)) ^«?/« w 

- A^--^K l/a+1 {(3r)K l/a {(3r)} , (51) 

1 roc ( fp 

(Tl{x)) Ren . = (Tlix))^. + — J o dppDl l=N (PR) l^K* /a+1 ([3r) 
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2£ + £) ^K l/a+l {(3r)K l/a {(3r) 



ar 



2 i K 2 l/a {(3r) 



(52) 



(Ti{x)) Ren . = (fi{x)) Reg . + ^j~ ' d(3(3D\ =N {(3R){{w 2 t 



P 2 



x 



a \ 4 



-K l/a+1 {(5r)K l/a {(3r) 



+ 



4^U-t + 



1\ , 4£ 7 (1 1 



ar 2 \a 2 



*?/«(/&•) , 



(53) 



and by boost invariance in the t — z plane {T^(x)) Ren , = (Tq (:r)) Ren .. Where 
the expression to (T^(x)) Reg , is given by 



( T ll( x )}Reg. = u 1 (a,7)diag(l, 1, -3, + - 1/6) u 2 (a, -f)diag(l, -1/2,3/2, 1) , 

(54) 

with 

UM ' 7) = 144Q7TV [(Q ' 4 " 1} " 30a " 47(1 " 7 ' )] (55) 



and 



u 2 (a,7) 



1 



2ttV 



i(a- 2 -l) + a- 2 7 ( 7 -i; 



(56) 



Substituting the asymptotic expression for D l n=N (f3R), Eq.s (31)- (35), into 
the above equations, after long calculation, we obtain: 



(t;{x)) 



Ren. 



(f»(x)) Reg . + 



1 



7 (R\ 2l/a (g? -i/af, 



N % 



X 



a 3 (a + 2 7 )(3a + 27) 



[(£ - l/6)dia^(ui,u 2 , «3,ui) 



7 



where 



+ o P - ^ ) diag(h 1 ,h 2 ,h 3 ,hi) 



vi = -(12a 3 + 4cry(8a + 77) + 87 s 
v 2 = Q a 3 + 2cry(5a + 27) , 



(57) 

(58) 
(59) 
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v 3 = - [6a 3 + 2a 7 (5a + 2 7 )] (3a + 2 7 ) , (60) 

and 

hi = 2a + 2 7 , (61) 

h 2 = 2a , (62) 

h 3 = -Q a - 4 7 . (63) 

From the above expressions we can see that: i) It reproduces the Guimaraes 
and Linet's result in the limit R — > 0, zi) the corrections decrease with the 
power of (R/r) proportional to the fractional part of </>/</>o and in) for £ = 1/6 
the trace of (57) vanishes. 

Finally it should be mentioned that this expression satisfies the conser- 
vation condition 

V v {T;{x))b,u. = . (64) 

As to the third model the respective Green function is given by (28) and 
the radial function a(r) by (11). Adopting a similar procedure as in the last 
two previous calculations we obtain: 

(T^x)) Ren . = (T°(x)) Reg . + —j o d(3(3 3 D ((3R){2^K 2 ((3r) 

+ 2^-^K^r)} , (65) 

I roo c 1 

{Tl{x)) Ren . = (Tl(x)) Reg . + — j o dW 3 D {m)[^Kl{Pr) 

- ^K^K^r) -Ix 2 (/3r)|, (66) 

(T*(x)) Ren . = (7^(x)) Bea . + — d(3[3 3 D ([3R) {2(£ - l/4)ff?(/?r) 

Air z a Jo 1 

+ ^K ((3r)K 1 ((3r) + 2^-l/A))K 2 ((3r)^ , (67) 
with {fl(x)) Ren , = (f$(x)) Ren .. Now we have 

(T;{x)) Reg . = y^tM' 4 - L-3,1) 

+ 7T7~5~t( Q; ~ 2 — 1)(^ — 1/Q)diag(2, —1, 3, 2) , (68) 



24ttV 



19 



which is the result for the renormalized vacuum expectation value of the 
energy-momentum tensor associated with a massless scalar field on a idealized 
cosmic string spacetime [14, 17]. This result, together with (39), seem very 
curious, because although we are considering a charged field in the presence 
of a cylindrical shell of magnetic field, differently from the two previous 
cases, the electromagnetic effect appears only on the seconds terms for both 
calculations; moreover this effect is of long-range type. 

Finally, substituting (42) into (65)-(67), taking into account the same 
considerations to overcome the integral problem on vacuum expectation value 
of the square of the scalar field, we obtain: 

(f;(x)) Ren . = (f;{x)) Reg . - A7r2ar ] Hqr) (Z - 1/6)^(2, -1, 3, 2) . (69) 

It is easy to see that this result is traceless when £ = 1/6 as conformal 
invariance demands on this spacetime; however as notice by Bruce Allen et 
al for £ 7^ 1/6, (69) presents a small violation of the conservation condition, 
which vanishes for R — > 0, or q — > oo. Being £ ^ 1/6, the energy density 
associated with a massless scalar charged field in this spacetime presents a 
long-range corrections arising from the magnetic field. For realistic models 
where the parameter a is very close to the unity this correction becomes even 
more relevant than the contribution coming from the geometry itself up to 
distance of the order the radius of the observable Universe. 



5. Concluding Remarks 

In this paper we have analyzed the vacuum polarization effects of a charged 
massless scalar field due to the locally flat cosmic string spacetime and a 
magnetic flux confined on a tube. 

Three distinct configurations of magnetic flux have been considered. Among 
them, the cylindrical shell of 5-function field presents a long-range effect on 
the vacuum polarization. 

In general, the results obtained to (T^(x)) Ren , as to (&*(x)&(x)) Ren ., out- 
side the tube, present corrections due to the magnetic tube which vanish 
when we take — > or R — > 0. 
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Specifically speaking about the trace anomaly, it is well known that [18] 

<t;mw = ^ , (to) 

where 02 (x) presents contributions coming from the geometry of the space- 
time as well as contribution from the gauge field strength [19]. However 
because the spacetime under consideration is locally flat, and the region of 
our interest has no magnetic field we have a,2(x) = 0. 

We conclude this paper by mentioning the importance to consider a finite 
structure of the magnetic flux on the quantum aspects of a charged field on 
an Abelian cosmic string spacetime. Although the complete knowledge about 
this structure could not be very well understood, its influence on the vacuum 
polarization effect can be even more relevant than the influence of geometry 
itself. 
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6. Figure Captions 



Figure 1: The solid curve in Fig. 1(a) represents the integrand of Eq. (30) 
written in terms of the new variable v = f3r for the first model considering the 
complete expression for the coefficient D\ given by (19). The dashed curve 
is the same integrand now using the approximated expression for Dq, Eq. 
(31). A similar analysis for the second model is given in Fig. 1(6). In both 
numerical analysis we have used a = 0.99, 7 = 0.2, N = and R/r = 10~ 3 . 

Figure 2: The solid curve in Fig. 2 represents the integrand of (39) for 
the third model considering the exact expression for D , given by (27). The 
dashed curve exhibits the same integrand considering now the approximated 
expression to D . In this numerical analysis we have used a = 0.99, 7 = 0.2, 
N = and R/r = 10~ 3 . 
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